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I describe recent developments in quenched chiral perturbation theory (QChPT) and the status of weak matrix 
elements involving light quarks. I illustrate how, with improved statistical errors, and with calculations of the 
masses of baryons containing non-degenerate quarks, there is now a clear need for extrapolations of higher order 
than linear in the quark mass. I describe how QChPT makes predictions for the functional forms to use in 
such extrapolations, and emphasize the distinction between contributions coming from chiral loops which are 
similar to those present in unquenched theories, and those from rj' loops which are pure quenched artifacts. I 
describe a fit to the baryon masses using the predictions of QChPT. I give a status report on the numerical 
evidence for rj' loops, concluding that they are likely present, and are characterized by a coupling S — 0.1 — 0.2. 
I use the difference between chiral loops in QCD and quenched QCD to estimate the quenching errors in a 
variety of quantities. I then turn to results for matrix elements, largely from quenched simulations. Results 
for quenched decay constants cannot yet be reliably extrapolated to the continuum limit. By contrast, new 
results for Bk suggest a continuum, quenched value of _B/f (NDR, 2 GeV) = 0.5977 ± 0.0064 ± 0.0166, based on a 
quadratic extrapolation in a. The theoretical basis for using a quadratic extrapolation has been confirmed. For 
the first time there is significant evidence that unquenching changes Bk, and my estimate for the value in QCD is 
-Bif (NDR, 2 GeV) = 0.66±0.02±0.11. Here the second error is a conservative estimate of the systematic error due 
to uncertainties in the effect of quenching. A less conservative viewpoint reduces 0.11 to 0.03. Non-perturbative 
renormalization has been successfully applied to four-fermion operators, and has been used to calculate Bk with 
Wilson fermions. The results are consistent with those for staggered fermions. Results for other matrix elements 
are mentioned briefly. 



1. INTRODUCTION 

This talk reviews the status of two distinct 
though related subjects: our understanding of 
chiral extrapolations, and results for weak ma- 
trix elements involving light (i.e. u, d and s) 
quarks. The major connection between these sub- 
jects is that understanding chiral extrapolations 
allows us to reduce, or at least estimate, the er- 
rors in matrix elements. Indeed, in a number of 
matrix elements, the dominant errors are those 
due to chiral extrapolation and to the use of the 
quenched approximation. I will argue that an 
understanding of chiral extrapolations gives us a 
handle on both of these errors. 

While understanding errors in detail is a sign 
of a maturing field, we are ultimately interested 
in the results for the matrix elements themselves. 
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The phenomenological implications of these re- 
sults were emphasized in previous reviews [^,||. 
Here I only note that it is important to calculate 
matrix elements for which we know the answer, 
e.g. /jr/Mp and Jk/ Itt, in order to convince our- 
selves, and others, that our predictions for un- 
known matrix elements are reliable. The reliabil- 
ity of a result for fjj, for example, will be gauged 
in part by how well we can calculate And it 
would be a real coup if we were able to show in 
detail that QCD indeed explains the A/ = 1/2 
rule in K ^ tttt decays. 

But to mc the most interesting part of the en- 
terprise is the possibility of using the lattice re- 
sults to calculate quantities which allow us to 
test the Standard Model. In this category, the 
light-quark matrix element with which we have 
had the most success is Bk- The lattice re- 
sult is already used by phenomenological anal- 
yses which attempt to determine the CP viola- 
tion in the CKM matrix from the experimental 
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number for e. I describe below the latest twists 
in the saga of the lattice result for Bk- What 
we would like to do is extend this success to the 
raft of B-parameters which are needed to predict 
e'/e. There has been some progress this year on 
the contributions from electromagnetic penguins, 
but we have made no headway towards calculat- 
ing strong penguins. I note parenthetically that 
another input into the prediction of e'/e is the 
strange quark mass. The recent work of Gupta 
and Bhattacharya and the Fermilab group |^ 
suggest a value of rric, considerably smaller than 
the accepted phenomenological estimates, which 
will substantially increase the prediction for e'/e. 

Much of the preceding could have been writ- 
ten in 1989, when I gave the review talk on weak 
matrix elements |^ . How has the field progressed 
since then? I see considerable advances in two 
areas. First, the entire field of weak matrix el- 
ements involving heavy-light hadrons has blos- 
somed. 1989 was early days in our calculation of 
the simplest such matrix elements, fjy and fs- In 
1996 a plethora of quantities are being calculated, 
and the subject deserves its own plenary talk |l0| . 
Second, while the subject of my talks then and 
now is similar, there has been enormous progress 
in understanding and reducing systematic errors. 
Thus, whereas in 1989 I noted the possibility of 
using chiral loops to estimate quenching errors, 
we now have a technology (quenched chiral per- 
turbation theory — QChPT) which allows us to 
make these estimates. We have learned that the 
quenched approximation (QQCD) is most prob- 
ably singular in the chiral limit, and there is a 
growing body of numerical evidence showing this, 
although the case is not closed. The reduction in 
statistical errors has allowed us to go beyond sim- 
ple linear extrapolations in light quark masses, 
and thus begin to test the predictions of QChPT. 
The increase in computer power has allowed us to 
study systematically the dependence of matrix el- 
ements on the lattice spacing, a. We have learned 
how to get more reliable estimates using lattice 
perturbation theory |lll|| . And, finally, we have 
begun to use non-perturbative matching of lat- 
tice and continuum operators, as discussed here 
by Rossi |l|. 

The body of this talk is divided into two parts. 



In the first. Sees. |[^, I focus on chiral extrapola- 
tions: why we need them, how we calculate their 
expected form, the evidence for chiral loops, and 
how we can use them to estimate quenching er- 
rors. In the second part. Sees. ^[^, I give an 
update on results for weak matrix elements. I 
discuss /^/Mp, /x/Ztt, Bk and a few related B- 
parameters. Results for structure functions have 
been reviewed here by Gockeler ||l^. There has 
been little progress on semi-leptonic form factors, 
nor on flavor singlet matrix elements and scatter- 
ing lengths, since last years talks by Simone jl^ 
and Okawa p^ . 



2. WHY DO WE NEED QChPT? 

Until recently linear chiral extrapolations (i.e. 
a + j3mq) have sufficed for most quantities. This 
is no longer true. This change has come about for 
two reasons. First, smaller statistical errors (and 
to some extent the use of a larger range of quark 
masses) have exposed the inadequacies of linear 
fits, as already stressed here by Gottheb ||l^. An 
example, taken from Ref. |0 is shown in Fig. |l|. 
The range of quark masses is ms/3 — 2tos, typi- 
cal of that in present calculations. While Ma is 
adequately fit by a straight line, there is definite, 
though small, curvature in Mjsi. The curves are 
the result of a fit using QChPT ||l|l, to which I 
return below. 

The second demonstration of the failure of lin- 
ear extrapolations has come from studying the 
octet baryon mass splittings |l^. The new fea- 
ture here is the consideration of baryons com- 
posed of non-degenerate quarks. I show the re- 
sults for (Ms — Mfi)/ [nis ~ rriu) (with m„ — rud) 
in Fig. Ij. If baryons masses were linear functions 
of TOs and ruu then the data would lie on a hor- 
izontal line. Instead the results vary by a factor 
of four. This is a glaring deviation from linear 
behavior, in contrast to the subtle effect shown 
in Fig. |. 

Once there is evidence of non-linearity, we need 
to know the appropriate functional form with 
which to extrapolate to the chiral limit. This is 
where (Q)ChPT comes in. In the second exam- 
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Figure 1. Mm and Ma versus M^, with quenched 
Wilson ferniions, at /3 = 6, on 32^ x 64 lattices. 
The vertical band is the range of estimates for 
rus. 
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Figure 2. Results for {M^ — Ma) / ims—mu) (dia- 
monds). The "burst" is the physical result using 
the linear fit shown. Crosses are from a global 
chiral fit. 



pie, the prediction is |19 1 

Ms - Ma 



SD ^ Ml - Ml „ M% - Ml 
6 rUs — rriu rUs — to„ 

+ei{mu + md + nig) + e'j(m„ + m^) , (1) 

where Mss is the mass of the quenched ss me- 
son.0 QChPT also implies that \d[ \ < \di\, which 
is why in Fig. || I have plotted the data ver- 
sus (M| — Ml)/{ms — niu)- The good news is 
that the points do lie approximately on a single 
curve — which would not be true for a poor choice 
of y-axis. The bad news is that the fit to the the 
di term and a constant is not that good. 

This example illustrates the benefits which can 
accrue if one knows the chiral expansion of the 
quantity under study. First, the data collapses 
onto a single curve, allowing an extrapolation to 
the physical quark masses. And, second, the the- 
oretical input reduces the length of the extrapola- 
tion. In Fig. ^, for example, to reach the physical 
point requires an extrapolation by less than a fac- 

have omitted terms which, though more singular in the 
chiral limit, are expected to be numerically small for the 
range of mq under study. 



tor of 2. This is much smaller than the ratio of 
the lightest quark mass to the physical value of 
(m„ -I- m£i)/2 — a factor of roughly 8. 

In fact, in the present example, the data are not 
yet accurate enough to distinguish between the d 
and e terms in Eq. |l|. A global fit of the QChPT 
prediction to this and other mass differences (the 
results of which are shown in the Figure) implies 
that both types of term are present . 

3. CHIRAL PERTURBATION THEORY 

This brings me to a summary of QChPT. In 
QCD, chiral perturbation theory predicts the 
form of the chiral expansion for quantities in- 
volving one or more light quarks. The expan- 
sion involves terms analytic in mq and the ex- 
ternal momenta, and non-analytic terms due to 
pion loops .0 The analytic terms are restricted in 
form, though not in magnitude, by chiral symme- 
try, while the non-analytic terms are completely 
predicted given the analytic terms. 

The same type of expansions can be developed 
in quenched QCD using QChPT. The method 
was worked out in Refs. EG-M, with the the- 



"Pion" here refers to any of the pseudo-Goldstono bosons. 
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oretically best motivated formulation being that 
of Bernard and Golterman [^. Their method 
gives a precise meaning to the quark-flow dia- 
grams which I use below. They have extended 
it also to partially quenched theories (those hav- 
ing both valence and sea quarks but with different 
masses) [Q. Results are available for pion prop- 
erties and the condensate |2^, Bk and related 
matrix elements Jb, Bb and the Isgur-Wise 
function |^,^ , baryon masses , and scat- 

tering lengths I will not describe technical 

details, but rather focus on the aims and major 
conclusions of the approach. For a more technical 
review see Ref. [|8| . 

As I see it, the major aims of QChPT are these. 

• To predict the form of the chiral expansion, 
which can then be used to fit and extrapolate the 
data. This was the approach taken above for the 
baryon masses. 

• To estimate the size of the quenching error 
by comparing the contribution of the pion loops 
in QCD and QQCD, using the coefficients of the 
chiral fits in QCD (from phenomenology) and in 
QQCD (from a fit to the lattice data). I return 
to such estimates in Sec. ^. 

I begin by describing the form of the chiral ex- 
pansions, and in particular how they are affected 
by quenching. The largest changes are to the 
non-analytic terms, i.e. those due to pion loops. 
There are two distinct effects. 

(1) Quenching removes loops which, at the under- 
lying quark level, require an internal loop. This 
is illustrated for baryon masses in Fig. |[ Dia- 
grams of types (a) and (b) contribute in QCD, 
but only type (b) occurs in QQCD. These loops 
give rise to terms in the chiral expansions. 
Thus for baryon masses, quenching only changes 
the coefficient of these terms. In other quantities, 
e.g. /,r, they are removed entirely. 

(2) Quenching introduces artifacts due to rj' 
loops — as in Fig. ^(c). These are chiral loops 
because the rj' remains light in QQCD. Their 
strength is determined by the size of the "hair- 
pin" vertex, and is parameterized by S and q;$ 
(defined in Sec. || below). 

The first effect of quenching is of greater prac- 
tical importance, and dominates most estimates 
of quenching errors. It does not change the form 
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(b) 
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Figure 3. Quark flow diagrams for Afj 
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of the chiral expansion, only the size of the terms. 

The second effect leads to new terms in the 
chiral expansion, some of which are singular in 
the chiral limit. If these terms are large, then 
one can be sure that quenching errors are large. 
One wants to work at large enough quark mass 
so that these new terms are numerically small. 
In practice this means keeping niq above about 
ms/4 - ms/3. 

I illustrate these general comments using the 
results of Labrenz and I for baryon masses . 
The form of the chiral expansion in QCD is^ 



(2) 



with C3 predicted in terms of g^NN and /^r. In 
QQCD 

M^..y = + ^2 Ml + c^Ml + c^M^+ 

6 (cf + c^Ml logM^) + a*c^Af3 + .(.3) 

The first line has the same form as in QCD, 
although the constants multiplying the analytic 
terms in the two theories are unrelated, is 
predicted to be non-vanishing, though different 
from C3. The second line is the contribution of 77' 

^There are also M^logM-^ terms in both QCD and 
QQCD, the coefficients of which have not yet been cal- 
culated in QQCD. For the limited range of quark masses 
used in simulations, I expect that these terms can be ad- 
equately represented by the terms, whose coefficients 
are unknown parameters. 
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Figure 4. Contributions to Af 



global chiral fit. All quantities in lattice units. 
Vertical lines indicate range of estimates of m^. 



loops and is a quenched artifact. Note that it is 
the dominant correction in the chiral limit. 

In order to test QChPT in more detail, I have 
attempted to fit the expressions outlined above to 
the octet and decuplet baryon masses from Ref. 
|T7[ . There are 48 masses, to be fit in terms of 
19 underlying parameters: the octet and decuplet 
masses in the chiral limit, 3 constants of the form 
c^, 6 of the form c^, 6 pion-nucleon couplings, 6 
and «$. I have found a reasonable description of 
the data with S « 0.1, but the errors are too large 
to pin down the values of all the constants |18[ . 
Examples of the fit are shown in Figs. |l| and 

Although this is a first, rather crude, attempt 
at such a fit, several important lessons emerge. 

(1) For present quark masses, one needs several 
terms in the chiral expansion to fit the data. This 
in turn requires that one have high statistics re- 
sults for a number of light quark masses. 

(2) One must check that the "fit" is not resulting 
from large cancellations between different orders. 
The situation for my fit is illustrated by Fig. |[ 
where I show the different contributions to M^. 
Note that the relative size of the terms is not de- 
termined by Mn, but rather by the fit as a whole. 
The most important point is that the terms 
are considerably smaller than those of 0{M^) up 



to at least niq = ttIs- The terms are part 
of a different series and need not be smaller than 
those of O(M^). Similarly the terms are the 
leading quenched artifact, and should not be com- 
pared to the other terms. Thus the convergence 
is acceptable for < TOs, though it is dubious 
for the highest mass point. 

(3) The artifacts (in particular the (5M^ terms) 
can lead to unusual behavior at small Af^r, as il- 
lustrated in the fit to Ma (Fig. p. 

(4) Since the "5-terms" are artifacts of quench- 
ing, and their relative contribution increases as 
Afjr ^ 0, it makes more sense phenomenologi- 
cally to extrapolate without including them. In 
other words, a better estimate of the unquenched 
value for A/a in the chiral limit can probably be 
obtained simply using a linear extrapolation in 
M^. This is, however, a complicated issue which 
needs more thought. 

(5) The output of the fit includes pion-nucleon 
couplings whose values should be compared to 
more direct determinations. 

(6) Finally, the fact that a fit can be found at 
all gives me confidence to stick my neck out and 
proceed with the estimates of quenching errors in 
baryon masses. It should be noted, however, that 
a fit involving only analytic terms, including up 
to Af^, can probably not be ruled out. 

What of quantities other than baryon masses? 
In Sec. 7.3 I discuss fits to Bk, another quantity 
in which chiral loops survive quenching. The data 
is consistent with the non-analytic term predicted 
by QChPT. Good data also exists for Mp. It 
shows curvature, but is consistent with either cu- 
bic or quartic terms |p9|. What do we expect from 
QChPT? In QCD the chiral expansion for Mp 
has the same form as for baryon masses . The 
QChPT theory calculation has not been done, but 
it is simple to see that form will be as for baryons, 
Eq. ||, except that = 0. Thus an M^ term is 
entirely a quenched artifact — and a potential win- 
dow on 

What of quantities involving pions, for which 
there is very good data? For the most part, 
quenching simply removes the non-analytic terms 
of QCD and replaces them with artifacts propor- 
tional to 5. The search for these is the subject of 
the next section. 
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Table 1 

Results from the quenched 77' two point function. 
W and S denote Wilson and staggered fermions. 



Rcf. 



Yr. 



a$ (3 W/S 



JLQCDp4| 
OSU| 
RomeTp6| 

osup 

FNAL 




0.14(01) 
0.27(10) 

« 0.15 
0.19(05) 

< 0.02 



0.7 5.7 W 
0.6 6.0 S 
5.7 W 
0.6 6.0 S 
>0 5.7 W 



4. EVIDENCE FOR if LOOPS 

The credibility of QChPT rests in part on 
the observation of the singularities predicted in 
the chiral limit. If such quenched artifacts are 
present, then we need to study them if only to 
know at what quark masses to work in order to 
avoid them! What follows in this section is an 
update of the 1994 review of Gupta Js^ . 

The most direct way of measuring 6 is from 
the T]' correlator. If the quarks are degenerate, 
then the part of the quenched chiral Lagrangian 
bilinear in the 77' is |23,Pql 



M'/- '/ (4) 
-(iV//3) {a^d^rj'd^W - mlrj'^) . (5) 



6 is defined by (5 = mg/(487r2/^)- In QQCD the 
terms in the second line must be treated as ver- 
tices, and cannot be iterated. They contribute to 
the disconnected part of the rj' correlator, whereas 
the first line determines the connected part. Thus 
the r]' is degenerate with the pion, but has ad- 
ditional vertices. To study this, various groups 
have looked at the ratio of the disconnected to 
connected parts, at p = 0, whose predicted form 
at long times is 



R{t) = t{Nf/3){ml - a^M^) /{2A'U) 



(6) 



I collect the results in Table |l|, converting mp into 
6 using a determined from nip. 

All groups except Ref. report a non-zero 
value for 6 in the range 0.1 — 0.3. What they 
actually measure, as illustrated in Fig. ^, is the 
combination ttiq— a$M^, which they then extrap- 
olate to = 0. I have extracted the results for 
a$ from such plots. As the figure shows, there is 
a considerable cancellation between the mo and 
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Figure 5. a^(mQ — aq,M^) from the OSU group. 



a$ terms at the largest quark masses, which cor- 
respond to A/tt ~ 0.8 GeV. This may explain why 
Ref. Q does not see a signal for i5. 

Clearly, further work is needed to sort out the 
differences between the various groups. As em- 
phasized by Thacker js^, this is mainly an issue 
of understanding systematic errors. In particu- 
lar, contamination from excited states leads to 
an apparent linear rise of R(t), and thus to an 
overestimate of mg. Indeed, the difference be- 
tween the OSU results last year and this is the 
use of smeared sources to reduce such contami- 
nation. This leads to a smaller 5, as shown in 
Fig. H Ref. also find that S decreases as the 
volume is increased. 

I want to mention also that the rf correlator 
has also been studied in partially quenched the- 
ories, with Nf = -6, -4, -2, and iV/ = 2,4 
p7[ |. The former work is part of the "bermion" 
program which aims to extrapolate from negative 
to positive Nf . For any non-zero Nf the analysis 
is different than in the quenched theory, because 
the hairpin vertices do iterate, and lead to a shift 
in the 77' mass. Indeed m^' is reduced (increased) 
for Nf < (> 0), and both changes are observed! 
This gives me more confidence in the results of 
these groups at Nf — 0. The bottom Hue appears 
to be that there is relatively little dependence of 
77^0 on Nf. 

Other ways of obtaining S rely on loop effects, 
such as that in Fig. |^(c). For quenched pion 
masses 77' loops lead to terms which are singular 
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in the chiral limit 



mi + m2 



+ -- 



1 - M log 



Mi, 



Mi, 



log 



Mi, 



Mti Mti 



A2 



+ C2(toi + m,) 



(7) 



Here Mij is the mass of pion composed of a quark 
of mass mi and antiquark of mass mj , A is an un- 
known scale, and c, an unknown constant. The 
tilde is relevant only to staggered fermions, and 
indicates that it is the mass of the flavor sin- 
glet pion, and not of the lattice pseudo-Goldstone 
pion, which appears. This is important because, 
at finite lattice spacing. Ma does not vanish in 
the chiral limit, so there is no true singularity. 

In his 1994 review, Gupta fit the world's data 
for staggered fermions at /3 = 6 having mi = m,- 
I have updated his plot, including new JLQCD 
data in Fig. ^ To set the scale, note that 
mgtt « 0.024. The dashed hue is Gupta's fit to 
Eq. ^, giving 5 — 0.085, while the solid line in- 
cludes also an m^ term, and gives 5 — 0.13. These 
non-zero values were driven by the results from 
Kim and Sinclair (KS), who use quark masses as 
low as O.lm^ |40|, but they are now supported by 
the JLQCD results. 

Last year, Mawhinney proposed an alternative 
explanation for the increase visible at small mq 
|4ll , namely an offset in the intercept of Af^ 



M(, = Co + n^{mi + m,) + 



(8) 



In his model, cq is proportional to the minimum 
eigenvalue of the Dirac operator, and thus falls as 
\/V. This model explains the detailed structure 
of his results for and (i/)-;/') at /3 = 5.7. It 
also describes the data of KS, except for volume 
dependence o/cq. As the Fig. ^ shows, the results 
of KS from 24^ and 32^ lattices are consistent, 
whereas in Mawhinney's model the rise at small 
mq should be reduced in amplitude by 0.4 on the 
larger lattice. 

The fit in Fig. ^ is for pions composed of de- 
generate quarks. One can further test QChPT by 



^The vertex leads to terms proportional to log M-^ 
which are not singular in the chiral limit, and can be rep- 
resented approximately by analytic terms. 
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Figure 6. Chiral fit to log(Af.^/m,) at /3 = 6. 
Some points have been offset for clarity. 

noting that Eq. ^ is not simply a function of the 
average quark mass — there is a predicted depen- 
dence on mi — m,. In Mawhinney's model, this 
dependence would presumably enter only through 
a (mi —m,)'^ term, and thus would be a weaker ef- 
fect. JLQCD have extensive data from the range 
P = 5.7 — 6.4, with both mi = m, and mi 7^ m,. 
They have fit to Eq. ^, and thus obtained S as 
a function of (3. They find reasonable fits, with 
S « 0.06 for most /3. I have several comments on 
their fits. First, they have used Ma rather than 
Mii, which leads to an underestimate of 5, par- 
ticularly at the smaller values of /3. Second, their 
results for the constants, particularly c,, vary 
rapidly with /3. One would expect that all dimen- 
sionless parameters in the fit (which are no less 
physical than, say, /^/Mp) should vary smoothly 
and slowly with p. This suggests to me that terms 
of O(m^) may be needed. Finally, it would be 
interesting to attempt a fit to the JLQCD data 
along the lines suggested by Mawhinney, but in- 
cluding a (mi — m,)'^ term. 

Clearly more work is needed to establish con- 
vincingly that there are chiral singularities in A/jr- 
One should keep in mind that the effects are 
small, ^ 5% at the lightest m,, so it is impressive 
that we can study them at all. Let me mention 
also some other complications. 
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(1) It will be hard to see the "singularities" with 
staggered fermions for /3 < 6. This is because 
Mii — Mii grows like (at fixed physical quark 
mass). Indeed, by /? = 5.7 the flavor singlet pion 
has a mass comparable to Mp\ Thus the 77' is 
no longer light, and its loop effects will be sup- 
pressed. In fact, the rise in M^/rriq as ^ 
for /3 = 5.7 is very gradual]!^, and could be due 
to the C2 term. 

(2) It will be hard to see the singularities using 
Wilson fermions. This is because we do not know, 
a priori^ where vanishes, and, as shown by 
Mawhinney, it is hard to distinguish the log di- 
vergence of QChPT from an offset in m^. 

(3) A related log divergence is predicted for (ipi/j), 
which has not been seen so far [^,^. It is not 
clear to me that this is a problem for QChPT, 
however, because it is difficult to extract the non- 
perturbative part of (f/'V') from the quadratically 
divergent perturbative background. 

Two other quantities give evidence concerning 
6. The first uses the ratio of decay constants 



Table 2 

Estimates of quenching errors. 



Rbg — fh/ifiih^) 



(9) 



This is designed to cancel the analytic terms 
proportional to niq p^ ], leaving a non-analytic 
term proportional to 6. The latest analysis finds 
S « 0.14 Q. It is noteworthy that a good fit 
once again requires the inclusion of 0{m'^) terms. 
The second quantity is the double difference 



ES2 = (Mn - Ma) - 3(AfH. - Ms-) 



(10) 



which is one measure of the breaking of the equal 
spacing rule for decuplets. This is a good window 
on artifacts due to quenching because its expan- 
sion begins at O(M^) in QCD, but contains terms 
proportional to SM^logA'U in QQCD The 
LANL group finds that ES2 differs from zero by 
2-cr ijl^, and I find that the data can be fit with 
<5«0.lll|. 

In my view, the preponderance of the evidence 
suggests a value of 6 in the range 0.1 — 0.2. All 
extractions are complicated by the fact that the 
effects proportional to S are small with present 
quark masses. To avoid them, one should use 
quark masses above ms/4 — ms/3. This is true 
not only for the light quark quantities discussed 
above, but also for heavy-light quantities such as 
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Ref. 



Error 



Jk/U - 1 
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0.2 
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-0.04 



0.05 
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0.3 



/b- This, too, is predicted to be singular as the 
light quark mass vanishes | p6[ |. 

5. QUENCHING ERRORS 

I close the first part of the talk by listing, in Ta- 
ble |2| a sampling of estimates of quenching errors, 
defined by 

I make the estimates by taking the numerator to 
be the difference between the pion loop contribu- 
tions in the full and quenched chiral expansions. 
To obtain numerical values I set A = nip (A is the 
scale occurring in chiral logs), use f = ^ Jk, 
and assume 6 — 0.1 and a$ = 0. For the es- 
timates of heavy-light quantities I set g' — 0, 
where g' is an rj'-B-B coupling defined in Ref. 
|p6|| . These estimates assume that the extrap- 
olation to the light quark mass is done linearly 
from TUq w mg/2. For example, fsa in QQCD 
is not the quenched value with the physical d- 
quark mass (which would contain a large artifact 
proportional to S), but rather the value obtained 
by linear extrapolation from ms/2, where the S 
terms are much smaller. This is an attempt to 
mimic what is actually done in numerical simula- 
tions. 

For the first two estimates, I have used the facts 
that, in QCD, [|| 

U « / [1 - 0.5L(Mk)] , (12) 
Ik/U « 1 - 0.25L{Mk) - 0.375L(M^) , (13) 
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(where L{M) = {M / An log(A'P / A"^ 
U = 93MeV), while in QQCD g|| 



Ik 



Ml 



Ml - Ml 



log 



M2 



Ml 



1 



and 

(14) 
(15) 



I have not included the difference of pion loop 
contributions to Mp, since the loop has not been 
evaluated in QChPT, and a model calculation 
suggests that the difference is small Details 
of the remaining estimates can be found in the 
references. 

Let me stress that these are estimates and 
not calculations. What they give is a sense of 
the effect of quenching on the contributions of 
"pion" clouds surrounding hadrons — these clouds 
are very different in QQCD and QCD! But this 
difference in clouds could be cancelled numeri- 
cally by differences in the analytic terms in the 
chiral expansion. As discussed in Ref. a 
more conservative view is thus to treat the esti- 
mates as rough upper bounds on the quenching 
error. Those involving ratios (e.g. fx/ f-n) are 
probably more reliable since some of the analytic 
terms do not contribute. One can also form dou- 
ble ratios for which the error estimates are yet 
more reliable (e.g. Rbg and ES2 from the previ- 
ous section; see also Ref. [^), but these quanti- 
ties are of less phenomenological interest. 

My aim in making these estimates is to obtain 
a sense of which quenched quantities are likely to 
be more reliable and which less, and to get an 
sense of the possible size of quenching errors. My 
conclusion is that the errors could be significant 
in a number of quantities, including those involv- 
ing heavy- light mesons. One might have hoped 
that the ratio fsJlB^ would have small quench- 
ing errors, but the chiral loops indicate otherwise. 
For some other quantities, such as BbJBb^ and 
Bk, the quenching errors are likely to be smaller. 

If these estimates work, then it will be worth- 
while extending them to other matrix elements of 
phenomenological interest, e.g. K — > tttt ampli- 
tudes. Then, when numerical results in QQCD 
are obtained, we have at least a rough estimate 
of the quenching error in hand. Do the estimates 
work? As we will see below, those for /ir/'Tip, 
fii/ f-n and Bk are consistent with the numerical 
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Figure 7. fi^/Mp with quenched Wilson fermions. 



results obtained to date. 

6. RESULTS FOR DECAY CONSTANTS 

For the remainder of the talk I will don the 
hat of a reviewer, and discuss the status of re- 
sults for weak matrix elements. All results will 
be quenched, unless otherwise noted. I begin with 
fir/Mp, the results for which are shown in Figs. 
(Wilson fermions) and || (SW fermions, with 
tadpole improved csw)- The normalization here 
is /^'^Pt = 0.131 MeV, whereas I use 93 MeV else- 
where in this talk. 

Consider the Wilson data first. One expects 
a linear dependence on a, and the two lines are 
linear extrapolations taken from Ref. |Q. The 
solid line is a fit to all the data, while the dashed 
curve excludes the point with largest a (which 
might lie outside the linear region). It appears 
that the quenched result is lower than experi- 
ment, but there is a 5 — 10% uncertainty. Im- 
proving the fermion action (Fig. |^) doesn't help 
much because of uncertainties in the normaliza- 
tion of the axial current. For the FNAL data, the 
upper (lower) points correspond to using as{n/a) 
{as{l/a)) in the matching factor. The two sets 
of UKQCD95 points correspond to different nor- 
malization schemes. Again the results appear to 
extrapolate to a point below experiment. 

It is disappointing that we have not done bet- 
ter with such a basic quantity. We need to reduce 
both statistical errors and normalization uncer- 
tainty. The latter may require non-perturbative 
methods, or the use of staggered fermions (where 
Za = 1). Note that chiral loops estimate that 
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the quenched result will undershoot by 12%, and 
this appears correct in sign, and not far off in 
magnitude. 

Results for {/k — J-n)/ f-n are shown in Fig. ||. 
This ratio measures the mass dependence of de- 
cay constants. Chiral loops suggest a 40% un- 
derestimate in QQCD. The line is a fit to all the 
Wilson data (including the largest a's), and in- 
deed gives a result about half of the experimental 
value. The new UKQCD results, using tadpole 
improved SW fermions, are, by contrast, rising 
towards the experimental value. It will take a 
substantial reduction in statistical errors to sort 
this out. 

7. STATUS OF Bk- STAGGERED 

is defined by 



Bk = 



{K\sjf^dsjf^d\K) 
{8/3){K\s^!:d\0){0\sjLd\K) 



(16) 



It is a scale dependent quantity, and I will quote 
results in the NDR (naive dimensional regulariza- 
tion) scheme at 2 GeV. It can be calculated with 
very small statistical errors, and has turned out 
to be a fount of knowledge about systematic er- 
rors. This is true for both staggered and Wilson 
fermions, though for different reasons. 



There has been considerable progress with both 
types of fermions in the last year. I begin with 
staggered fermions, which hold the advantage 
for Bk as they have a remnant chiral symme- 
try. Back in 1989, I thought we knew what the 
quenched answer was, based on calculations at 
/3 = 6 on 16^ and 24^ lattices: Bk = 0.70(2) 
1^,^. I also argued that quenching errors were 
likely small (see Table ||). I was wrong on the 
former, though maybe not on the latter. 

By 1993, Gupta, Kilcup and I had found 
that Bk had a considerable a dependence Q. 
Applying Symanzik's improvement program, I 
argued that the discretization errors in Bk 
should be O(a^), and not 0{a). Based on 
this, we extrapolated our data quadratically, and 
quoted BK{NDR,2GeV) = 0.616(20)(27) for 
the quenched result. Our data alone, however, 
was not good enough to distinguish linear and 
quadratic dependences. 

Last year, JLQCD presented results from a 
more extensive study (using /3 = 5.85, 5.93, 6 
and 6.2) |5^. Their data strongly favored a lin- 
ear dependence on a. If correct, this would lead 
to a value of Bk close to 0.5. The only hope 
for someone convinced of an dependence was 
competition between a number of terms. 

Faced with this contradiction between numer- 
ical data and theory, JLQCD have done further 
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work on both fronts They have added two 

additfonal lattice spacings, (3 = 5.7 and 6.4, thus 
increasing the lever arm. They have also carried 
out finite volume studies at /? = 6 and 6.4, find- 
ing only a small effect. Their data are shown 
in Fig. |l^. "Invariant" and "Landau" refer to 
two possible discretizations of the operators — the 
staggered fermion operators are spread out over 
a 2"* hypercube, and one can either make them 
gauge invariant by including gauge links, or by 
fixing to Landau gauge and omitting the links. 
The solid (dashed) lines show quadratic (linear) 
fits to the first five points. The x^/d.o.f. are 



Fit Invariant 
a 0.86 



a 



1.80 



Landau 
0.67 
2.21 



thus favoring the linear fit, but by a much smaller 
difference than last year. If one uses only the 
first four points then linear and quadratic fits are 
equally good. What has changed since last year 
is that the new point at /3 = 6.4 lies above the 
straight line passing through the next four points. 

JLQCD have also checked the theoretical argu- 
ment using a simpler method of operator enumer- 
ation!^,^ The conclusion is that there cannot 
be 0{a) corrections to Bk, because there are no 
operators available with which one could remove 

® A similar metliod has also been introduced by Luo |57 . 



these corrections. Thus JLQCD use quadratic ex- 
trapolation and quote (for degenerate quarks) 

Bif(NDR, 2 GeV) = 0.5977±0.0064±0.0166 , (17) 

where the first error is statistical, the second due 
to truncation of perturbation theory. This new 
result agrees with that from 1993 (indeed, the 
results are consistent at each /3), but has much 
smaller errors. To give an indication of how far 
things have come, compare our 1993 result at f3 — 
6 with Landau-gauge operators, 0.723(87) Q, to 
the corresponding JLQCD result 0.714(12). 

The perturbative error in Bk arises from trun- 
cating the matching of lattice and continuum op- 
erators to one-loop order. The use of two different 
lattice operators allows one to estimate this er- 
ror without resort to guesswork about the higher 
order terms in the perturbative expansion. The 
difference between the results from the two oper- 
ators is of 0[a(2 GeV)^], and thus should remain 
finite in the continuum limit. This is what is ob- 
served in Fig. |l^. 

I will take Eq. |l^ as the best estimate of Bk 
in QQCD. The errors are so much smaller than 
those in previous staggered results and in the re- 
sults with Wilson fermions discussed below, that 
the global average is not significantly different 
from the JLQCD number alone. The saga is not 
quite over, however, since one should confirm the 

dependence by running at even smaller lattice 
spacings. JLQCD intend to run at /3 = 6.6. If 
the present extrapolation holds up, then it shows 
how one must beware of keeping only a single 
term when extrapolating in a. 

7.1. Unquenching Bk 

To obtain a result for QCD proper, two steps 
remain: the inclusion of dynamical quarks, and 
the use of nis ^ rud- The OSU group has made 
important progress on the first step ||3l|] . Previous 
studies (summarized in Ref. |^) found that Bis- 
was reduced slightly by sea quarks, although the 
effect was not statistically significant. The OSU 
study, by contrast, finds a statistically significant 
increase in Bk 



BK(NBR,2GeY,Nf = 3) 
BK(NDR,2GeV, AT/ = 0) 



= 1.05 ±0.02. 



(18) 
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They have improved upon previous work by re- 
ducing statistical errors, and by choosing their 
point at Nf = {(3 = 6.05) to better match the 
lattice spacing a.t Nf — 2 {f3 — 5.7, mqU — 0.01) 
and 4 (/? = 5.4, m^a = 0.01). 

There are systematic errors in this result which 
have yet to be estimated. First, the dynami- 



^sca 



T val 



and so they are truely unquenched sim- 
ulations. But mP*^^^ is determined by extrapolat- 
ing in the valence quark mass alone, and is thus 
a partially quenched result. This introduces an 
uncertainty in rus which feeds into the estimate 
of the Nf dependence of Bk ■ Similarly, a is de- 
termined by a partially quenched extrapolation, 
resulting in an uncertainty in the matching fac- 
tors between lattice and continuum operators. 

But probably the most important error comes 
from the possibility of significant a dependence 
in the ratio in Eq. |l^. The result quoted is for 
= 2GeV, at which a the discretization error 
in the quenched Bk is 15%. It is not inconceiv- 
able that, say, Bk in QCD has very little depen- 
dence on a, in which case the ratio would increase 
to ~ 1.2 in the continuum limit. Clearly it is very 
important to repeat the comparison at a different 
lattice spacing. 

Despite these uncertainties, I will take the OSU 
result and error as the best estimate of the effect 
of quenching at a = 0. I am being less conserva- 
tive than I might be because a small quenching 
error in Bk is consistent with the expectations 
of QChPT. A more conservative estimate for the 
ratio would be 1.05 ± 0.15. 

7.2. Bk for non-degenerate quarks 

What remains is to extrapolate from to^ 
nid ~ to rris — m^^'^^ and nid 

This appears difficult because it requires dynam- 
ical quarks with very small masses. This may not 
be necessary, however, if one uses ChPT to guide 
the extrapolation |Q. The point is that the chi- 
ral expansion in QCD is p8| , p3[ 

Bk 
B 



phys 

- my . 



= 1- 3 



y\ny + by + cye^ 



(19) 



where 

e = (m, - md)/(m, +md), y = M|./(47r/)2,(20) 



and B, b and c are unknown constants. At this 
order / can be equally well taken to be or /k- 
Equation |l^ is an expansion in y, but is valid for 
all e. The idea is to determine c by working at 
small e, and then use the formula to extrapolate 
to e = 1. This ignores corrections of 0{y^), and 
so the errors in the extrapolation are likely to be 
- 25%. 

Notice that m„ does not enter into Eq. |l^. 
Thus one can get away with a simulation using 
only two different dynamical quark masses, e.g. 
setting ruu — ra^ < m^^^^ /2, while holding nis + 
TUd = m^^y^. To date, no such calculation has 
been done. To make an estimate I use the chiral 
log alone, i.e. set c = 0, yielding 

BK(non-degen) = (1.04 - 1.08)SK(degen) .(21) 

The range comes from using f = f-n and Jk, 
and varying the scale in the logarithm from nip — 
1 GeV. Since the chiral log comes mainly from 
kaon and rj loops |23|, I prefer / = /a', which 
leads to 1.04 — 1.05 for the ratio. To be conser- 
vative I take 1.05 ± 0.05, and assume that the 
generous error is large enough to include also the 
error in the estimate of the effect of unquenching. 
This leads to a final estimate of 

5a(NDR, 2 GeV, QCD) = 0.66±0.02±0.03 , (22) 

where the first error is that in the quenched value, 
the second that in the estimate of unquenching 
and using non-degenerate quarks. Taking the 
more conservative estimate of the unquenching 
error (15%), and adding it in quadrature with 
the (5%) estimate of the error in accounting for 
non-degenerate quarks, increases the second error 
in Eq. || to 0.11. 

It is customary to quote a result for the renor- 
malization group invariant quantity 



B 



K 



1 



Itt 



/3i7o - /3o7i 



in the notation of Ref. [|^. Using as(2GeV) = 
0.3 and iV/ = 3, I find Bk = 0.90(3)(4), with 
the last error increasing to 0.14 with the more 
conservative error. This differs from the result I 
quoted in Ref. Q, because I am here using the 
2-loop formula and a continuum choice of a^. 
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7.3. Chiral behavior of Bk 

Since Bk can be calculated very accurately, 
it provides a potential testing ground for (par- 
tially) quenched ChPT. This year, for the first 
time, such tests have been undertaken, with re- 
sults from OSU JLQCD H, and Lee and 
Klomfass |59|. I note only some highlights. 

It turns out that, for e = 0, Eq. |l^ is valid 
for all Nf This is why my estimate of the 

quenching error for Bk with degenerate quarks in 
Table | is zero. Thus the first test of (P)QChPT 
is to see whether the —Sylny term is present. 
The OSU group has the most extensive data as 
a function of y, and indeed observe curvature of 
the expected sign and magnitude for Nf = 0, 2, 4. 
JLQCD also finds reasonable fits to the chiral 
form, as long as they allow a substantial depen- 
dence of / on lattice spacing. They also study 
other B parameters, with similar conclusions. 

Not everything works. JLQCD finds that the 
volume dependence predicted by the chiral log 
|23| is too small to fit their data. Fitting to the 
expected form for e ^ in QQCD, they find 
S = —0.3(3), i.e. of the opposite sign to the 
other determinations discussed in Sec. ^. Lee 
and Klomfass have studied the e dependence with 
Nf = 2 (for which there is as yet no PQChPT pre- 
diction). It will be interesting to see how things 
evolve. My only comment is that one may need 
to include 0{y^) terms in the chiral fits. 

8. STATUS OF Bk- WILSON 

There has also been considerable progress in 
the last year in the calculation of Bk using Wil- 
son and SW fermions. The challenge here is to 
account for the effects of the explicit chiral sym- 
metry breaking in the fermion action. Success 
with Bk would give one confidence to attempt 
more complicated calculations. 

The operator of interest, 

Ov+A = 's'^^^ds-ff^d + 57/^75^^ 57^75^ , (23) 
can "mix" with four other dimension 6 operators 

OHa = Zv+A (^v+A + ^ z,o}j + 0(a) (24) 
where the O on the r.h.s. are lattice operators. 



The Oi are listed in Refs. j6^j6^. The meaning of 
this equation is that, for the appropriate choices 
of Zv+A aud the Zi, the lattice and continuum 
operators will have the same matrix elements, up 
to corrections of 0{a). In particular, while the 
matrix elements of a general four fermion opera- 
tor has the chiral expansion 

{K\0\K)^a + (3Ml + 5iMli+ (25) 

PK-pKh + S2M]^ + 53PK-pK) + ■■■ , (26) 

chiral symmetry implies that a — P = 6i = for 
the particular operator O = Oy"^. Thus, one 
can test that the Zi are correct by checking that 
the first three terms are absent.^ Note that the 
Zi must be known quite accurately because the 
terms we are removing are higher order in the 
chiral expansion than the terms we are keeping. 

Five methods have been used to determine the 
Zi and Zy+A- 

(1) One-loop perturbation theory. This fails to 
give the correct chiral behavior, even when tad- 
pole improved. 

(2) Use (1) plus enforce chiral behavior by ad- 
justing subsets of the Zi by hand j63|. Different 
subsets give differing results, introducing an ad- 
ditional systematic error. Results for a variety of 
a were presented by Soni last year |^ . 

(3) Use (1) and discard the a, /3 and 7 terms, 
determined by doing the calculation at a variety 
of momenta. New results using this method come 
from the LANL group ||6^ . Since the Zi are incor- 
rect, there is, however, an error of 0{g^) in Bk- 

(4) Non-perturbative matching by imposing con- 
tinuum normalization conditions on Landau- 
gauge quark matrix elements. This approach has 
been pioneered by the Rome group, and is re- 
viewed here by Rossi The original calcula- 
tion omitted one operator ]6C| ] , but has now been 
corrected [ |62[ . 

(5) Determine the Zi non-perturbatively by im- 
posing chiral ward identities on quark matrix el- 
ements. Determine Zy+A as in (4). This method 
has been introduced by JLQCD . 

The methods of choice are clearly (4) and (5), 
as long as they can determine the Zi accurately 

have ignored chiral logarithms, which will complicate 
the analysis, but can probably be ignored given present 
errors and ranges of Mj^. 
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Figure 11. Quenched Bk with Wilson fermions. 

enough. In fact, both methods work well: the 
errors in the non-perturbative results are much 
smaller than their difference from the one-loop 
perturbative values. And both methods find that 
the matrix element of Oy"^\ has the correct chiral 
behavior, within statistical errors. What remains 
to be studied is the uncertainty introduced by 
the fact that there are Gribov copies in Landau 
gauge. Prior experience suggests that this will be 
a small effect. 

It is not yet clear which, if either, of methods 
(4) and (5) is preferable for determining the Zi. 
As stressed in Ref. |^ the Zi are unique, up to 
corrections of 0{a). In this sense, both methods 
must give the same results. But they are quite 
different in detail, and it may be that the errors 
are smaller with one method or the other. It will 
be interesting to see a comprehensive comparison 
between them and also with perturbation theory. 

In Fig. |ll| I collect the results for Bk- All cal- 
culations use TOs = iTid and the quenched approx- 
imation. The fact that most of the results agree is 
a significant success, given the variety of methods 
employed. It is hard to judge which method gives 
the smallest errors, because each group uses dif- 
ferent ensembles and lattice sizes, and estimates 
systematic errors differently. The errors are larger 
than with staggered fermions mostly because of 
the errors in the z,-. 



Extrapolating to a = using the data in Fig. 
[n] would give a result with a large uncertainty. 
Fortunately, JLQCD has found a more accurate 
approach. Instead of Bk, they consider the ra- 
tio of the matrix element of Oy^\ to its vac- 
uum saturation approximant. The latter differs 
from the denominator of Bk (Eq. |6|) at finite 
lattice spacing. The advantage of this choice is 
that the Zi appear in both the numerator and 
denominator, leading to smaller statistical er- 
rors. The disadvantage is that the new ratio has 
the wrong chiral behavior at finite a. It turns 
out that there is an overall gain, and from their 
calculations at (3 = 5.9, 6.1 and 6.3 they find 
BA:(NDR,2GeV) = 0.63(8). This is the result 
shown at a = in Fig. ^ It agrees with the 
staggered result, although it has much larger er- 
rors. Nevertheless, it is an important consistency 
check, and is close to ruling out the use of a linear 
extrapolation in a with staggered fermions. 

9. OTHER MATRIX ELEMENTS 

The LANL group Q has quenched results (at 
(3 — 6) for the matrix elements which determine 
the dominant part of the electromagnetic penguin 
contribution to e'/e 



B 



1=3/2 _ 



B 



1=3/2 



= 0.58±0.02(stat) 
= 0.81 ± 0.03(stat) 



+0.07 
-0.03 

+om 

-0.02 



(27) 
(28) 



These are in the NDR scheme at 2 GeV. The sec- 
ond error is from the truncation of perturbation 
matching factors. These numbers lie at or below 
the lower end of the range used by phenomenolo- 
gists. The LANL group also finds Bd = 0.78(1). 

There are also new results for fp and [^|9) , 
for the pion polarizability |66|, and for strange 
quark contributions to magnetic moments ||67[| . 

10. FUTURE DIRECTIONS 

This year has seen the first detailed tests of 
the predicted chiral behavior of quenched quan- 
tities. Further work along these lines will help 
us make better extrapolations, and improve our 
understanding of quenching errors. It is also a 
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warm-up exercise for the use of chiral perturba- 
tfon theory in unquenched theories. I have out- 
hned one such appUcation in Sec. 7.2. I expect 



the technique to be of wide utiUty given the dif- 
ficulty in simulating light dynamical fermions. 

As for matrix elements, there has been substan- 
tial progress on Bk- It appears that we finally 
know the quenched result, thanks largely to the 
efforts of JLQCD. At the same time, it is disturb- 
ing that the complicated a dependence has made 
it so difficult to remove the last 20% of the er- 
rors. One wonders whether similar complications 
lie lurking beneath the relatively large errors in 
other matrix elements. 

The improved results for Bk with Wilson 
fermions show that non-perturbative normaliza- 
tion of operators is viable. My hope is that we 
can now return to an issue set aside in 1989: 
the calculation oi K ^ tttt amplitudes. The 
main stumbling block is the need to subtract 
lower-dimension operators. A method exists for 
staggered fermions, but the errors have so far 
swamped the signal. With Wilson fermions, one 
needs a non-perturbative method, and the hope 
is that using quark matrix elements in Landau 
gauge will do the job j62|. Work is underway 
with both types of fermion. Given the success of 
the Schrodinger functional method at calculating 
current renormalizations ||6^, it should be tried 
also for four fermion operators. 

Back in 1989, 1 also described preliminary work 
on non-leptonic D decays, e.g. D Ktt. Almost 
no progress has been made since then, largely be- 
cause we have been lacking a good model of the 
decay amplitude for Euclidean momenta. A re- 
cent proposal by Ciuchini et al. may fill this gap 
|69|. 

Enormous computational resources have been 
used to calculate matrix elements in (P)QQCD. 
To proceed to QCD at anything other than a 
snail's pace may well require the use of improved 
actions. Indeed, the large discretization errors in 
quenched staggered Bk already cry out for im- 
provement. The fact that we know Bk very ac- 
curately will provide an excellent benchmark for 
such calculations. Working at smaller values of 
the cut-off, 1/a, alleviates some problems while 
making others worse. Subtraction of lower di- 



mension operators becomes simpler, but the eval- 
uation of mixing with operators of the same di- 
mension becomes more difficult. It will be very 
interesting to see how things develop. 
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